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We provide a theoretical framework to compute the zero field muon spin relaxation rate of a 
Heisenberg ferromagnet at low temperature. We use the linear spin wave approximation. The rate, 
which is a measure of the spin lattice relaxation induced by the magnetic fluctuations along the easy 
axis, allows to estimate the magnon stiffness constant. 
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I. INTRODUCTION 



Recently, the measured temperature behavior of the zero field muon,s.ein relaxation (//SR^-fate of some Heisenberg 
ferromagnets has been analysed quantitatively in the paramagnetictTu and ferromagnetica critical regimes. This 
analysis has given information on the anisotropy of the spin dynamics in the reciprocal and direct spaces. It has shown 
that the anisotropy in the reciprocal space, which is induced by the dipolar interaction, has a strong influence on the 
\Q • dynamics near the zone center of the Brillouin zone. Note that this interaction always exists in real ferromagnets. The 
Q\ ' experimental data in_the isotropic dipolar Heisenberg regime are well explained by the mode coupling approximation 
of Frey and Schwabl.El On the other hand, up to now, the low temperature behavior of the relaxation rate has only been 
considered qualitativelyfl or in restricted physical cases.Lj The main purpose of this work is to provide a theoretical 
framework to analyse ^SR relaxation data recorded on Heisenberg ferromagnets at low temperature. 

The organization of this paper is as follows. In Sec. EJ, we summarize the r elati on between the measured depolar- 



in 

ization function and the spin-spin correlation function of the magnet. In Sec. Ill we compute the muon spin lattice 
relaxation rate in the linear spin wave approximation for a Heisenberg hamiltonian with an energy gap. In Sec. IV 
we apply the result of the previous section to describe the spin lattice relaxation of a dipolar Heisenberg magnet in 
I . the spin wave temperature region. In this last section we present the conclusion of our work. 
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II. MUON SPIN RELAXATION AND SPIN-SPIN CORRELATION FUNCTIONS 

O ■ 

We take the Z-axis parallel to the incoming muon beam polarization. A zero field muon spin relaxation measurement 
, consists of measajxiag the time depolarization function Pz (t) .Em For simplicity we take the easy magnetic axis z parallel 
to Z. We writeEM 

a- 

Pz(t) = exp[-V*(i)], (la) 

with 

/•* 

„2 



Mt) =% dr{t~ r) [<Z> xx (t) + * w (r)] . (lb) 



7 M is the muon gyromagnetic ratio : 7^ = 851.6 Mrad s~ 1 T~ 1 . i/j z (t) does not depend on the muon pulsation frequency 
w M because the associated energy Uuj^ is negligible In Eq. (|fy) <& qq (t) is the symmetrized time correlation function 
of the a fluctuating component of the local magnetic field at the muon site and {x, y, z} is an orthogonal frame. 

Experimentally it is usually observed that Pz(t) is an exponential function. Eq. ([!]) predicts such a result if r is 
very small relative to t, i.e. if the characteristic decay time of $ QQ is much smaller than t which is ~ 2 (is. Taking 
this hypothesis which is justified and the fact that the field correlation functions are even functions of r, we derive 
that Pz{t) is an exponential function characterized by a damping rate X z which can be written in terms of the time 
Fourier transform of field correlation functions at u) = 

K =Tzil[® xx (w = Q) + % y (w = 0)]. (2) 
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It is possible to express the magnetic field at the muon site as a function of, on one hand, a tensor which describes 
the coupling between the localized spins of the metal and the muon spin and on the other hand of the localized spins 
components themselves. Because 3> qq (t) is quadratic in fields, A z can be written as a linear combination of spin-spin 
correlation functions of the magnet. In general the calculation of these correlation functions is performed in q space, 
i.e. in the first Brillouin zone. Therefore the spatial Fourier transform of these functions has to be introduced. This 
is done in detail in Ref. [l| Using the result given at Eq. (g) we find 



A, 



|^£[G^(q)G-(- 



-q) 



+G^(q)G a7 (~q)] A^(q) 



(3) 



The integral extends over the first Brillouin zone. This equation shows that A 2 depends on the coupling between 
the muon spin and the spins of the magnet through the coupling tensor G(q) and the spin correlation tensor of the 
magnet, A(q, w), taken at u> = 0, i.e. A(q) = A(q, uj — 0). The information on the muon localization site is contained 
in G(q). We have defined V = (h /4tt) 7^ (<?lA<b) where fj,o is the permeability of free space, <7l the Lande factor 
and hb the Bohr magneton. V is the volume of the sample. The formula given at Eq. (||) allows to deal with Bravais 
and non Bravais lattices, for instance the hexagonal close packed lattice. In the latter case G(q) = l/id^^tq) 
where the index d runs over the rid non equivalent sites of the magnet and Gd(q) = exp[iq • (i + d)]G ri+d . i runs 
over each cell of the crystal lattice and G Ti+d is a dimensionless tensor which accounts for the classical dipolar and 
Fermi contact couplings between the muon spin and the magnet spins located at distance vector r^+d from the muon. 
A(q) writes J^d d' ^-dd'(q) where Add'(q) is the correlation tensor between spins belonging to cell sites d and d' . Note 
that the case of Bravais lattice is included in the definition we have just given : d takes then only one value (rid = 1) 
and the sums reduce to a single element. 

The tensor G(q) is the sum of dipolar and hyperfine tensors which we note D(q) and H(q) respectively. In general 
terms, the sublattice d contribution to D(q) writeso 



4tt P^(q)-G7(q) 



a/9 1 



(4a) 



where P^(q) is the longitudinal projection operator and C^ p (q) a symmetric tensor : 



a/3. 



(q) = 



1 — exp 



4q 



K^O 



+ 2(^)3/2 [ 2 ^ r i+d,o t n+d,l3(P3/2 (Af+d) ~ (ff^f+d)] ex P " r *+d) ■ 



(4b) 



The -D(q) expression is derived using an Ewald transformation.EHI The ip m (x) functions are defined in Ref. |l|. K is a 
vector of the reciprocal lattice and (respectively ro+d) is the vector that links the muon site to the ion belonging 
to sublattice d and located in lattice cell i (respectively origin) of the crystal lattice. Expression (4b) gives the same 
result for all values of the Ewald parameter g, but for numerical applications a value of g is chosen which ensures that 
both series ofEq. @ converge rapidly. The Gd(q) tensor, the trace of which is 1, reveals the symmetry of the point 
group at the muon site. Whereas the elements of the Gd(q) tensor are analytical functions of q for all q values, the 
elements of Pl(<1) are only piecewise continous at q = 0. 

The hyperfine interaction is short range and usually isotropic. In the lowest order in q we have 



H a P(q = 0) = l/n d Y,Hf{ci = 0) = r^HS* 



(5) 



where r M is the number of nearest neighbor magnetic ions to the muon site and H a constant which can be deduced 
from the muon spin rotation frequency at low temperature. Equation (|5|) is derived under the hypothesis that the 
muon site is a center of symmetry. 

In order to proceed further, in this paper we will only consider the behavior of G(q) near the zone center. From 
the previous results, using G Q ^(q) = l/^Ed^'d' 3 ^)' we derive in the limit q — > 



G ali {q) = -4tt 



G a/3 (0) 



4tt 



(6) 
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III. MUON SPIN RELAXATION RATE IN AN HEISENBERG FERROMAGNET WITH A SMALL 

ENERGY GAP 



For the computation of the damping rate we need an expression for the correlation functions of the magnet. 
We first have to specify its hamiltonian. As we have done beforeQ we suppose that the magnet is described by a 
Heisenberg interaction with a small energy gap. Because the minimum magnon energy is much larger than Hlo^, 
the energy conservation principle tells us that only the parallel (to the easy axis z) fluctuations contributes to the 
depolarization,El i.e. the measurements probe only the correlation function A zz (q). Then a close look at Eq. (||) 
indicates that the diagonal terms of G(q) do not contribute to the depolarization. This means that an isotropic 
hyperfine interaction does not influence the measured damping rate. For many possible muon localization sites in 
crystals we have C zx (q = 0) = C zy (q = 0) = 0. Using this hypothesis we deduce the simple result G z @(q — > 0) = 
-AnPf{q) for = x, y. This leads to the important fact that, within our hypothesis, the muon spin relaxation rate 
is independent of the muon localization site. The previous analysis leads to the following simple result : 



K = — dVf A"(q). (7) 
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A generalization of this result when C zx (q — 0) and C zy (q = 0) are not zero is not difficult : it results in a modifr 
prefactor that depends on C zx (q = 0) and C zy (q = 0). We compute A zz (q) in the linear spin wave approximation. 
In this approximation the fluctuating part of the z component of the total angular momentum of the magnetic ion, 
<5Jq, is expressed in terms of a sum of products of boson operators a J and a^: 

K = 4^ i (" q+k - ki ' K ^ 1 (8) 

k,ki 



where K is a vector of the reciprocal lattice. Each of the k and ki sums extends over N vectors of the first Brillouin 
zone. Because <5Jq is given in terms of products of two operators, the computation of A zz {q) requires the evaluation 
of four operator products. We decouple these products. We set : 

< a k a k! a q a qi > = Ar 2 <5 kiq <5 kqi 7l k (n q + 1), (9) 

where n q is the standard Bose occupation factor. With these results it is easy to derive an expression for the parallel 
spin-spin correlation function. We find : 

In order to proceed further we need an expression for the magnon dispersion relation hujq. The simplest choice is to 
write 

frujq = D m q 2 + A, (11) 

where D rn is the magnon stiffness constant and A is the energy gap of the magnon dispersion relation. We have 
supposed that the dispersion relation is isotropic at small q. The sum in Eq. ( |l0| ) can be converted to an integral 
and can be computed analytically using Eq. ([fl]). We obtain : 

A zz (q) VH ^ 



16^ 2 Dig 



exp[p m g 2 /4 + A) /k B T] - 1 
1 



exp [(D m q 2 BZ + A) /k B T] - 1 



(12) 



where k B is the Boltzmann constant and q B z the radius of the Brillouin zone. We note that usually at low temperature 
we have (D m q BZ + A) 3> k B T. This is the case, for instance, for GdNis for which D m — 3.2 meV.A 2 and q B z = 0.9 
A" 1 (see Ref. |). Then the expression of the correlation function greatly simplifies : 



A«(g) = VK 

{q> 16^ Dlq 



X 



expKA^ 2 /4 + A)/k B T]-l 



(13) 
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Notice that, because of the simple relation we take for the magnon dispersion relation, the correlation function is a 
function of the magnitude of q and not of its orientation. 

So far in this section we have not mentioned if we are dealing with Bravais or non Bravais lattices. In a non Bravais 
lattice where A zz (q) writes J^d d , A dd , (q), the existence of more than one atom per unit cell gives rise to acoustic and 
optical branches in the dispersion relation. Because of the relatively huge gap present in the optical branch(es), the 
spin correlation functions k z dd ,{q) with d ^ d! are not relevant as seen at Eqs. ( |l0| ) or (p^). 

The expression of the muon spin relaxation rate for the model considered in this section (specified by the dispersion 
relation given at Eq. (pi])) is obtained using Eqs. (@) and @. We derive 



D 3 



P(q BZ ), 



(14a) 



with 



P(q) = In 



1 - exp [-(£>„ 



<7 2 /4 + A)/k B T] 



1 - exp (-A/k B T) 



(14b) 



where we have defined C 
A) > D m q 2 BZ /A > k B T, we have 



(2/IStt) (/xo^tOSX^I 



129.39 (meV) 3 .A 6 .s _1 .K- 2 . In the limiting case, {D m q 2 BZ /A- 



A, 



CgjT 2 
Dl 



In 



exp (A/fcgT) 
exp (A/k B T) - 1 



(15) 



A numerical study of the P{q) function shows that the expression for \ z given at Eq. ( |l5| ) is a reasonable approximation 
of the A z expression given at Eq. ([l4]) . 

In practice we have A <C k B T. In this high temperature regime we obtain the simple result : 



A 2 



CgjT 2 



\n{k B T/A). 



(16) 



The T 2 \n(k B T / A) dependence of Aj has been predicted in the past for an anisotropic contact interaction which can 
occur in nuclear magnetic resonance. oEj Some years ago we have derived this temperature dependence for muon spin 
relaxation rate using the model considered in the present work.@ But the derivation was not sound mathematically : 
the constant G (introduced at Eq. (15) of Ref. ( |q ) was undetermined. 

The physical origin of the T 2 factor in Eq. (Hq) is clear : each of the two magnons contributing to the muon spin 
depolarization process (Raman process) accounts for a factor T. We note that A z given by Eq. (^6|) diverges if A = 
0, %.£-. for a pure Heisenberg magnet. This fact is not disturbing because a real magnet has always a small energy 
gap £3 

Although the result of Eq. ( |l6| ) has not been published before, it has already been used to extract the stiffness 
constant of the dipolar ferromagnet GdNis-from its zero muon spin relaxation rate dataO 

X z has been computed by other authorsQ for EuO. In the present work, we have presented a general framework to 
compute the /iSR relaxation rate for a ferromagnet at low temperature. Then for a given magnon dispersion relation 
the rate has been computed explicitely. In the proposed framework the coupling between the muon spin and the 
electronic spins is described by the tensor G(q). For a given muon localization site in a crystal, its elements are easily 
computed. One does not have to consider separately the relaxation due to the hyperfine and dipolar interaction. 



IV. MUON SPIN RELAXATION IN A DIPOLAR HEISENBERG FERROMAGNET AT LOW 

TEMPERATURE 

In the previous section we have computed A z for a simple Heisenberg hamiltonian with an energy gap. In a dipolar 
magnet the gap is induced by the dipolar interaction between the ions. Recently the neutron scattering function 
has been computed for a dipolar Heisenberg magnet in the linear spin wave approximation.!^ From this work an 
expression for the zz correlation function can be deduced. Relative to our expression the differences are that a 
relatively complicated weighting factor exists besides the thermal factor in Eq. (|Io| ) and the magnon dispersion 
relation for a dipolar Heisenberg magnet is used. However, qualitatively this factor does not influence the correlation 
function A zz (q). The magnon dispersion relation writeaiij 
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^1 = eq(eq + esin 2 9q), 



(17) 



where # q is the polar angle of wave vector q relative to the easy axis, £ q the dispersion relation without the dipolar 
interaction and e = (u, g L fi B M )/?i a characteristic dipolar energy (M is the saturation magnetization). Therefore 
a natural choice is to identify the gap A with he. In this case the effect of sin 2 # q is neglected. An other possibility is 
to take for the gap value the bulk magnetic anisotropy energy : A = gLHBB a where B a is the magnetic anisotropy 
field. These two choices represent limiting cases. We note that for GdNi 5 we have Ue/ks ~ 1-2 K and (gLl^BB a )/kB 
~ 0.3 K, i.e. there is a factor 4 difference. But because the purpose of the measurement of the relaxation rate is to 
measure the spin wave stiffness constant which is proportionnal to [ln^^T/A)] 1 / 3 , this uncertainty by a factor of 4 
will only introduce an error of ~ 20 % on the extracted D m value from the relaxation data (in the case of GdNis the 
D m value is deduced from the relaxation data recorded at about 10 K - 15 K). 

In this work we have given a general framework to analyse positive muon longitudinal relaxation data at low tem- 
perature for a Heisenberg magnet. Using the simplest possible approximation, we have shown that the measurements 
allow to estimate the magnon stiffness constant. The framework which is given here should allow to deal with other 
problems in an easy way. We think for instance of antiferromagnets. 
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